Anisotropic Node Removal in d-wave Superconductors under Magnetic Field 
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A phenomenological model that considers different secondary id xy gap amplitudes and quasipar- 
ticle effective charges for each nodal direction, is proposed to explain the observed anisotropic node 
removal by a magnetic field in high-T c cuprates. Two independent parity-breaking condensates 
develop, implying the induction of a magnetic moment per each nodal direction. The 
model outcomes are in agreement with the experimentally found relation Aj ~ \f~B. The secondary 
gap vanishes through a second-order phase transition at a critical temperature whose value for 
underdoped nodal-oriented YBCO is estimated to be ~ 6K for a field of 15T. 

PACS numbers: 74.20.Rp; 74.20.-z; 74.72.-h 
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Thanks to a wealth of experiments and theoretical 
studies on high T c cuprates, it has become clear that any 
viable theory of high-T c superconductivity should incor- 
porate, beside the superconducting and normal phases, 
several other ordered phases that may be externally in- 
duced or spontaneously generated, and can either coexist 
or compete. 

The idea of the existence of competing ground states 
and quantum critical points in high-T c cuprates was used 
by Votja, Zhang and Sachdev 0,0 to explain an observed 
anomalously large damping rate of quasiparticles near 
the gap nodes of the BiiSriC 'aC^Og+s superconductor. 
The quasiparticle damping rate, studied in photoemis- 
sion spectroscopy (ARPES) and optical conductivity 
experiments, revealed a linear in temperature depen- 
dence, in contrast to the cubic dependence expected in 
the conventional d-wave picture [5j . This anomalous en- 
hancement of the rate was considered in Ref. 1, ,2j to be 
associated to a quantum critical behavior and interpreted 
as the result of the fluctuations of a new, secondary order 
parameter that removed the quasiparticle nodes. Then, a 
phenomenological theory of nodal quasiparticles coupled 
to the fluctuations of a secondary (is or id xy ) pairing or- 
der parameter was proposed, and evidence of a quantum 
critical point in the superconducting phase was provided 
within a perturbative e-expansion in Ref. . This quan- 
tum critical point was later established beyond pertur- 
bation theory by Khveshchenko and Passke 6] , finding a 
damping rate ~ T at high temperatures. 

The development of a new order parameter has been 
also found in tunneling experiments with Y ' BaiCu-sO-i - x 
(YBCO) 0. In these studies, a splitting of the zero bias 
conductance peak (ZBCP) is observed, when a magnetic 
field is perpendicularly applied to the CuC>2 planes of 
the YBCO films. A ZBCP is associated to the presence 
of nodes, as it reflects the existence of Andreev surface 
bound states at and near the Fermi surface [8j. Hence, 
the splitting of the ZBCP may be an indication of the 
appearance of a secondary gap, and therefore of node re- 
moval by the magnetic field. In addition to the field in- 
duced splitting of the ZBCP, a spontaneous -probably as- 
sociated to the doping level- splitting was also seen in . 



These findings motivated Dagan and Deustscher |£| to 
perform a comprehensive study of the doping dependence 
of the ZBCP splitting in YBCO samples with and with- 
out external magnetic field. These authors concluded 
that above certain doping level near optimum doping, 
a spontaneous ZBCP splitting occurs, which increases 
with doping. It also increases with magnetic field applied 
along the c-axis. In underdoped samples the magnetic 
field could induce the ZBCP splitting, but no sponta- 
neous splitting was observed. 

Although the ZBCP splitting by the magnetic field 
may be considered as due to node removal, a different ex- 
planation, based on Doppler shift of the Andreev bound 
states has been also offered [lOj. However, a new set of 
experiments 

.11] 

has definitely favored node removal in 
the case of underdoped materials, measured in decreasing 
fields to eliminate hysteresis effects. In these experiments 
a large ZBCP splitting was clearly observed in films with 
thickness smaller than the London penetration depth. 
That is, under conditions where the Doppler shift effect is 
negligible. Another important finding of Ref. 0] i s that 
the splitting is very anisotropic. Only for samples having 
(1,1,0) orientation, the splitting was substantial. The au- 
thors of 0] noticed that the experimental data obtained 
for the (1,1,0) oriented samples were in good agreement 
with the predictions made by Laughlin |l2| in a model he 
had proposed some years ago to explain the behavior of 
the thermal conductivity of cuprates in a magnetic field 
[l3j . However, they also called attention that at present 
there is no theoretical model that predicts the observed 
anisotropy. Laughlin model, which is based on a transi- 
tion from a pure d-wave order parameter to a parity and 
time-reversal symmetry violating d x 2_ y 2 + id xy ground 
state, removes the nodes of both nodal directions. 

It should be mentioned that the behavior of the ther- 
mal conductivity of cuprates with a magnetic field has 
proven itself to be anisotropic too. The hallmark of Kr- 
ishana et al. effect |l3|, namely, the decreasing of the 
thermal conductivity with the magnetic field, followed 
by a kink and a plateau region, is indeed observed in 
only one of the two nodal orientations ^ij. This strik- 
ing result was revealed, after a long debate, by careful 
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experiments realized by Ando et al. [14|. Anisotropic 
in-plane transport properties of cuprates have been also 
observed by Ando and coworkers in measurements of the 
electrical resistivity of underdoped samples and of 
the magnetoresistance of lightly doped ones [16j . These 
experimental findings have indicated that the anisotropy 
between the nodal directions is related to the presence 
of charge stripes along the nodal directions of these ma- 
terials. Thus, the electrons organize in an anisotropic 
macroscopic state, characterized by stripes that depend 
on the temperature and doping. 

In this paper we propose a phenomenological model for 
the low-energy quasiparticle excitations at the nodes of a 
superconductor on which both magnetic field and doping 
can induce the transition from the regular d-wave pair- 
ing to d x 2_ y 2 + id xy - or d x 2_ y 2 + is-pairings. The novel 
feature of the present model is that it incorporates the 
observed anisotropy of the nodal directions, so it can de- 
scribe the anisotropic field-induced ZBCP splitting phe- 
nomenon 0, llfij . Two main new ideas have been imple- 
mented here: (1) to consider the induction of a secondary 
order parameter with id xy (or is) symmetry that may 
have different amplitudes along the two nodal directions, 
and, (2) inspired by the above mentioned anisotropic 
charge distribution along the nodal directions, to allow 
for different effective charges of the quasiparticle excita- 
tions about these directions. 

Our treatment will follow the standard proce- 
dure based on the rotationally covariant Nambu 
representation, employed in Refs. 0, 0] and @], to 
describe the quasiparticles in terms of two species 
of Nambu bispinors. Then, the bispinors are com- 
bined in four-component Dirac spinors, which permit 
to write the action in a familiar field theory form. 
Following Lee and Wen [l7| . we assume that the 
spectrum of the quasiparticle excitations in the super- 
conductor is given by E(k) = ye(k) + A 2 (k). Here 
e(k) = 2t(cos(bk x )+cos(bk y )), as usual. For the proposed 
complex order parameter (c^c_^) = A(k), however, 
we propose a secondary id xy component characterized 
by two amplitudes Ai xy , (i — 1,2 ) , one for each nodal 
direction, given by A(k) = Ao(cos(6fca;) — cos(bk y )) + 



iAi xy sm(bk x )sin(bk y )[sin 2 (^(k x - ki x ))sin 2 (^(k y 
k ly ) + sin 2 (^(k x + k lx ))sin 2 {\{k y + k ly ))} + (1 <-> 2), 
where ±fej = ±(kp,±kp) are the momentum of the 
nodes along the two nodal directions. A similar order 
parameter can be used to produce is secondary pairing, 
by substituting the factor $,vn(bk x )sin(bk y ) by 1 in 
A(A:). From now on we concentrate our discussion on 
the id xy , as the results can be easily extended to the 
is case. Our guidance to propose such a form of the 
secondary order parameter has been to use a basis 
function that transforms under the symmetry group of 
the superconductor C iv x Z 2 18] in the same way as the 
usual id xy basis function. 

Rotating the coordinate system by 7r/4 and expanding 



about each nodal point up to linear order, one can write 
the low energy effective action of the quasiparticle exci- 
tations around the nodes as the sum of a fermion kinetic 
term 



S v = i ^ / d3;z #ia (lod + v F j x d x + v A j y d v ) tfia 
a =] J 

+ (lf la -> (fi2a,X < ► y)],(l) 

a Higgs- Yukawa (HY) term 

S v< p = J d 3 x[g<t) 1 Tp la ipi a + gfo ^2 V2<x&<x], (2) 



and a scalar field action, whose generic form can be 
found, after integrating out high-energy fermion modes, 
to be 
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In the above equations m, A, g and c are undetermined 
parameters to be fixed by the experiment, vf — 2y/2tb 
and va — are the characteristics velocities of the 
d-wave superconductor, and gfa = Ai xy and g(f>2 = 
A 2 xy We are using the irreducible representation 7 M = 
(t2, iri, — its), fi = 0,1,2, of the (2+l)-dimensional 
Dirac algebra, where Ti are the Pauli matrices. The 
quasiparticle states near the pairs of the opposite nodes 
have been described in terms of the bispinors (pi a = 
(c Q (hi) ,s af 3cl {-hi)), with lp 2a = ^^p^a, and we in- 
troduced the conjugates Tp la = (p| Q 7o, <? 2q = ^2q7o- The 
index a =\, \ denotes the two spin components, and e a jj 
is an antisymmetric tensor with en = 1. 

To make the fermion sector of the action explicitly 
Lorentz invariant, we can make the variable change x — > 
y, in the terms with the ipi a field, and 

— y, in those with the field (f2a- In 




this way, v — ^/vfva plays the role of the speed of light. 
From now on, for simplicity, we set all the velocities to 
1. 

At this point it is convenient to consider a reducible 
4x4 representation of the Dirac algebra by introducing 
the matrices T M = diagfap, — 7/*)) and combining the 
bispinors of each nodal direction in four-component Dirac 
spinors, to form two four-spinors of different "flavors" 
i = 1,2 defined as \&i = {<px\-, { p\\)\ ^2 = {j-f*i\i ¥*2j.) with 
conjugate \ti = ^To. Notice that the i th Dirac spinor 
involves the fields of the i th nodal direction only. In the 
4x4 representation the action becomes 
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where A = ir rir 2 = diag(h,-h)- The action Q is 
invariant under parity transformations 



(a; -x,y -> y) , <^ 



(5) 



with Tp = Ti (£> n . 

Let us now switch on a magnetic field along the c-axis. 
As mentioned above, an important point of the proposed 
model is to consider that the effective charges of the two 
nodal directions are different. Thus, the magnetic field 
couples with different strengths e\ and e 2 to fields and 
^>2 respectively. Then, in the presence of the magnetic 
field the action Q takes the form 

r 2 _ _ 

S B = if^^r^+^fiAtj 
J i=i 

+~(^i) 2 -™ 2 ^-A§] (6) 

where Z?f = <9 M — ieiA 1 ^*, and we take a gauge on which 
Aff* = (0, 0, Sxi) . If due to the magnetic field the scalar 

field acquires nonzero vacuum expectation value (vev) fa 
(induced secondary energy gap), the corresponding i th - 
parity symmetry is broken. To explore such a possibility 
it is needed to minimize the free energy of the system 



with respect to fa. I n the lowest approximation, it gets 
contributions from all tadpole diagrams of the effective 
theory. At finite temperature the minimum equation is 
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where (3 1 = T, and, to simplify the notation, we 
dropped the use of the "hat" in the scalar vev's. 

Since this is an effective low-energy theory of quasipar- 
ticle excitations around the nodes, the low-energy domain 
will be limited by the induced gap. Therefore, the diver- 
gent integrals should have cutoff Aj. Moreover, it is nat- 
ural to expect a magnetically induced gap smaller than 
the field generating it. As a consequence, the fermion 
infrared dynamics will be dominated by the lowest Lan- 
dau level, and the leading contribution to the fermion 
tadpole, second term in the RHS of (J7J), will come from 
1 = 0. Thereby, after doing the temperature sums, we 
obtain 
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In JSJ we took m 2 = 0, and introduced the notation 

a = "J^-- Eq. iJS} can have different field-dependent 
solutions, for different ranges of A and g. However, to 
describe the experimental findings, we need to consider 
the parameter range that leads to Aj ~ \/\eiB\ fllj |. 
Such a solution can be found neglecting the first term in 
the LHS of iJSJl, which corresponds to neglect the contri- 
bution of the tree-level Xfa term in the free-energy. In 
this way, the zero-temperature minimum results 



Aj 0) = 2^72^51 



(9) 



Experiments with YBCO samples gives the ratio -4= = 

i y B 

l.lmeV /T 1 / 2 in data taken in decreasing fields Re- 
instating the superconductor characteristic velocity in 
(jSJ through the change etB — ► (hv 2 /c)eiB, using Hv ~ 
0.25eVA for YBCO, and = e for one of the nodal di- 
rections we find a ~ 0.13. It can be checked that this 
result is not in contradiction with the dropping of the 



cubic term in (JSJ. 

We underline that the solution is the only extremal 

at zero temperature. That is, A^ = is not even an 
extremal of the free-energy at T = 0. Moreover, the 
breaking of the symmetry by the magnetic field at T = 
occurs independently of the strength of the Yukawa 
coupling <7, i.e. there is no critical g in this case. 

A similar zero temperature behavior was previously 
found in our studies of a (3+l)-dimensional HY model in 
the presence of a magnetic field 0,113] • However, there 
is a fundamental difference between that case and the one 
we are considering here. In the model of Refs. flfl licj. 
the magnetic field breaks a discrete chiral symmetry pro- 
ducing a chiral condensate (^^). Such a symmetry 
breaking is a manifestation of the phenomenon, intro- 
duced by Gusynin, Miransky, and Shovkovy [ij, of mag- 
netic catalysis of chiral symmetry breaking (MCCSB) 
[22| . The MCCSB mechanism has been suggested as a 
possible explanation for several transport properties of 
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strongly correlated electronic systems 23] . Nevertheless, 
in the id xy theory considered in the present paper, the 
condensation phenomenon has a different physical char- 
acter, since the Lagrangian has explicitly broken chi- 
ral symmetry due to the HY term. As it happens in 
MCCSB, here too the fermion infrared dynamics in the 
presence of a magnetic field induces the appearance of 
a gap and a fermion- antifcrmion condensate, but in this 
case it is a parity-breaking condensate instead 
of the chiral condensate ('J'S) associated to the MCCSB. 
We call this new mechanism of magnetic-field-induced 
dynamical symmetry breaking: magnetic catalysis of par- 
ity symmetry breaking. Physically, the parity condensate 
corresponds to the induction of a magnetic moment by 
the magnetic field. The induction of a magnetic moment 
was proposed by Laughlin as the key to explain the 
odd behavior of the thermal conductivity at B =^ in 
high-T c superconductors 0,0. 

As it is known, thermal effects tend to erase the gen- 
erated condensates. In the present case, this is re- 
flected in the tendency of the gap solution to vanishing 
(Aj(T) — > 0) as T is increased. That is, we find that at 
some critical temperature T C i, that depends on the mag- 
netic field, a second order phase transition takes place 
from a gapped to a gapless state for the i th nodal di- 
rection, independently of what occurs in the other nodal 
direction. Near the transition temperature T C i the gap 
solution can be written as 

A i (T) = ^(T? i -T 2 ) (10) 

with T cl = a^/\e l B\ = { v ^/8)A ( i 0) . Using the gap 

Aj- -* = 4meV", which corresponds to a field of 15T in the 
600 A film of undcrdoped (l,l,0)-oriented YBCO, mea- 
sured in decreasing field yjj , we estimate a critical tem- 
perature T c ~ 6K. 

Finally, we would like to comment on the anisotropic 
behavior of the proposed model in the presence of a mag- 
netic field. As seen from Eq. ©, the induced gap Ai(T) 
and the critical temperature T C i of the i th nodal direction 
depend on the strength of the effective charge of the 
quasiparticles of that particular direction. The smaller 
the effective charge e^, the smaller the Aj(T) and the 
T c j. Then, given a magnetic field B, we have that for 
T c2 {B) < T < T cl {B) no ZBCP splitting will be ob- 
served if the film is oriented along direction 2, but it 
will be present for films oriented along 1. In this sense, 
it would be interesting to have available data of ZBCP 
measurements taken alon g th e nodal direction perpen- 
dicular to the one used in [llj, as in that case one could 
use the relation e 2 = (^) 2 ^i to estimate the ratio of 
the two effective charges. We think that this anisotropic 
condensation mechanism is also the basis to explain the 
anisotropic transport properties in underdoped supercon- 



ducting samples [H H E . 
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